Integrability and ideal conductance at finite temperatures 
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_v Abstract 

a 
o 

O ■ We analyse the finite temperature charge stiffness D{T > 0), by a gener- 

•y-i . alization of Kohn's method, for the problem of a particle interacting with a 

X 

3 . fermionic bath in one dimension. We present analytical evidence, using the 

Bethe ansatz method, that D[T > 0) is finite in the integrable case where the 

mass of the particle equals the mass of the fermions and numerical evidence 

that it vanishes in the nonintegrable one of unequal masses. We conjecture 

that a finite D{T > 0) is a generic property of integrable systems. 
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This work relates the finite temperature charge transport to the response of energy levels 
to an infinitesimal fiux and proposes a connection to the integrability of quantum systems. 

Starting with the formulation by Kohn [|I|, the charge stiffness (or Drude weight) repre- 
senting the weight of the delta function contribution D5{ijj) part in the dynamical conduc- 
tivity o"(co'), has been investigated as a criterion for a metallic, superconducting or insulating 
state 1|,§]. The original approach [^ relates D to the response of the ground state energy to 
a magnetic fiux and thus requires only the calculation of an equilibrium property, bypassing 
a complete evaluation of the Kubo formula @] . 

A similar concept appears also in the study of the conductance of disordered metallic 
systems. Starting with the work of Thouless a close relation has also being established 
between the conductance, the sensitivity on boundary conditions (being equivalent to the 
introduction of a fiux) and the statistical properties of single-particle energy spectra of 
metallic systems. 

Recently it has also been observed that the level statistics in many-particle correlated 
systems is closely related to the integrability of the system |0,|^. Still the implications of 
this fact on transport quantities in correlated fermion systems at finite temperatures T > 0, 
e.g. a nonzero d.c. resistivity p{T) > or on the contrary a possible finite charge stiffness 
D(T) > 0, have not been known so far. 

In this work we analyse the relation between integrability and transport in the generic 
problem of a single tagged particle moving in a bath of fermions by a generalization of 
Kohn's method at finite temperatures. 

Kohn's method atT>0: We consider a general tight-binding Hamiltonian of the form 

H = -tJ2{e''^4_,,d, + H.c.) + H,nt = f + H,nt, (1) 

i 

representing a ID system of length L with periodic boundary conditions pierced by a fiux 
L(f), using the Peierls construction. T is the kinetic and Hmt the interaction part of the 
Hamiltonian. From the Kubo formula ||^,0] we can relate the imaginary part of the dynamical 
conductivity cr"{uj) to the charge stiffness D, 



D = -[ua {uj)]^^o = jiT^i-T) - 2^ Pn —— }, (2) 

where j is the current operator, (T) the thermal expectation value of the kinetic energy and 
Pn = e~^"/"^/Z the Boltzmann weight for an eigenstate | n) of Hamiltonian (||) with energy 

en- 

On the other hand |1]| we can evaluate, by second order perturbation theory, for (p —>■ 
a shift of the level | n) : 

^ ^ , (77 7 777) 1 -N 

e„(0) = (n I ff(0 = 0) I n) - <l){n \ j \ n) - <{>' J2 \ _^ - ^^^^^ I ^ I ^)- (3) 

Extracting second order terms in (f) (the curvature of levels) we see that: 

1 1 ld'^en{<p) ,,. 

d^-y:p„d„^-i:p„-^^. (4) 

In this demonstration we assumed that the levels | n) are nondegenerate but the same 
procedure can be used for degenerate levels after applying first order degenerate perturbation 
theory to lift the level degeneracy. This expression for D{T) generalizes Kohn's method 
to T > and reduces its evaluation to the calculation of curvature of levels under an 
external flux. Easier than the complete evaluation of Kubo formula, it can be performed 
either numericaly on finite size systems followed by finite size scaling or in some cases (e.g. 
integrable systems) analytically |]10|- It also provides a basis for analysing the conditions for 



the occurence of ideal conductance. 

The ground state of the metal should be characterized hy Dq = D(T = 0) > ||l]-0|, 
where Dq is a measure of of the charge carriers coherent motion. On physical grounds it 
is commonly assumed that D(T > 0) = 0, for generic macroscopic interacting systems 
(involving Umklapp scattering), except in the superconducting phase. This statement, never 
really proved to our knowledge for interacting fermions, implies according to our Eqs.(0,^ 
that the positive curvature of the ground state Do is cancelled by (mostly) negative curva- 
tures of low lying states, a situation reminiscent of the mechanism for phase transitions and 
the dissapearance of long range order. 



Model system: The Hamiltonian describing a particle in a bath of spinless fermions is: 
H = -thY.i^"^4+A + H.c.) -tY,{cUiC, + H.c.) + UY.dU4c,, (5) 

i i i 

where Ci{cl) are annihilation (creation) operators for A^ spinless fermions and di{dl) for the 
(tagged) particle on an L site chain with periodic boundary conditions. The interaction 
conies only through the on-site repulsion U > 0. The current j in (^ refers to the particle 
only ( j = —ithYjidi+idi + H.c) and [j, -ff] 7^ 0. The volume normalization 1/L in (^) is 
absent as the current refers only to one particle. 

The model (|^) is well adapted to our study: it is integrable by the Bethe ansatz method 
in the case of equal masses, i.e. th = t [|1^] (equivalent to the problem of a Hubbard chain 
in a nearly polarized state S^ = S^^^ — 1) and nonintegrable for unequal masses (th 7^ t). 
We have also previously studied Dq [O] of the particle as well as its quasiparticle properties 



UK. 

We first present the results for D{T) from exact numerical diagonalization of the Hamil- 
tonian (^ on finite size systems. We will consider only the A^ = L/2 case, which in this 
system should be analogous to all other fillings. By numerical calculation of the energy 
spectrum with and without a small flux (typically ~ 10~^) we deduce the energy level 
curvature from: 



1 en(0) + en(- 0) - 26(0) 
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In Fig.l we present D{T) as a function of th/t, for different temperatures T, size systems 
(from 6 to 14 sites) but fixed U/t = 2. We observe: a) a nonmonotonic behavior of D{T) as 
we sweep through the integrable th = t point for T > 0, b) a very weak dependence of D(T) 
with system size for the integrable case th = t while a rather strong one for the nonintegrable 
ones. 

These observations are the first hints that something particular happens at the integrable 
point. To further study this behavior, we present in Fig. 2 D(T)/Do as a function of 1/L for 
th = t. We indeed observe a scaling in 1/L which we also know that it is the relevant one 



to recover the ground state Dq |jT2|. The circle at 1/L = indicates the result of the Bethe 
ansatz calculation for th = t in very good agreement with the extapolated numerical results. 

The same plot in Fig.2 of In D{T)/ Do as a function of L^ for th = 0.5, 1.2, 2.0t seems to 
indicate a strong L dependence following a law D{T) ~ e~'^^ for these nonintegrable cases. 
This is consistent with a rapid decrease of D{T) for systems larger than the mean free path. 

Next to obtain an impression of the overall behavior of D{T) as function of temperature 
we present in Fig.3 D{T)/Dq for th = t,0.5t and different size systems. For th = 0.5t 
D(T)/Do seems to scale to zero for any T > as we found above while for th = t there 
is a smooth decrease with T. At high temperatures D{T) oc 1/T, vanishing for T — > cxd a 
general feature of systems on a lattice. Finally the discrepancy from the free particle D(T) 
indicates that the temperature effects cannot be accounted for by a simple scaling of Dq. 

Bethe ansatz analysis: We present here an analytical approach for the calculation of 



D{T) for the integrable case using the Bethe ansatz method along the line of reference [|T2 



The starting point is the Lieb- Wu [0 solution of the Hubbard model adapted for one spin up 
fermion (the tagged particle) and A^ spin down fermions. The Bethe ansatz wavef unctions, 
in the presence of flux 0, are then characterized by A^ + 1 quantum numbers kj given by the 
following equations (we take th = t = 1): 

Lkj = 27iIj + 9{sinkj-A), j = 1,...,N + 1 (7) 

e{p) = -2tan~\4p/U) (8) 

7V+1 7V+1 

Lj2kj = 2nJ2 h + 27r J + L0. (9) 

Every state is characterized by a set of half-odd integers {/j} and the (half-odd) integer 
J for (even) odd number of fermions. However this set of states(regular) does not constitute 
a complete set; to obtain a complete set we must include states representing bound states of 
energy of order U. The simplest way to include these states is by considering an electron-hole 
transformation (q = cj, di = (— l)*(ij) and then solving Eq.(|^ but for a L — N number of 
fermions (equal to the number of holes). 



The total energy of a regular state is given by : 



N+l 

E= Ee(^,) = E(-2cosfc,), (10) 



and for a bound state: 



E = U+ Y. ^2 COS kj. (11) 

i=i 

In the following we present explicitely only the analysis for the contribution of regular 
states as it is clear how to include the bound ones by the prescription described above |T^ 
(although we include both contributions in the results presented in Fig. 2 and 3). 

To order 1/L and for UL ^ 1 the values of kj are given by: 

k^=k'^ + Le{smk'^-A), A;° = ^. (12) 

To order one then the energy and momentum of a state can be written as: 

E = Y, <k^j) + 7 sin fc° ^(sin A;° - A) (13) 

iE^(sinA;°-A) = ^ + 0. (14) 

In this scheme the interaction between the particle and the fermions is represented 
through second term in Eq.(0), a correlation energy; the coupling to the flux is through the 
collective coordinate A. In the thermodynamic limit, defining a density p{k) (— tt < k < +tt) 
we obtain: 

E{p{k),A) = — f dkp{k){~2cosk + ^smke{smk-A)) (15) 

2tt J L 

^ [ dkp{k)9{smk-A) = P + (l), P=^. (16) 

To calculate the thermal average of the curvature of energy levels we will work in the 
grand canonical ensemble for the system plus the one particle. Assuming that the equilibrium 
distribution of k's is not affected by the presence of the one particle we use the Fermi-Dirac 
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distribution for free fermions f{k) = 1/(1 + exp((/i — t{k))/T) (yU is the chemical potential). 
We then assume that the distribution of A's is determined by the average correlation energy 

ec(A), 

e^(A) = — / dkf{k)2 sin k ^(sin A; - A), (17) 

through a Boltzmann weight: w{K) = exp(— ec(A)/T). 
Finally D{T) is given by: 



D = ^ J^^ dAg{A)w{A)^ J dkfik)D{A,k), (18) 



with g{A) = dP/dA determined from Eg. ([161): 

Z^ = ^jdAg{A)w{A), (20) 

(in the total partition function Z\ we add the contribution from both the regular and bound 
states). 

D(A,k) = -2s.nki^(-)'' + -^). (21) 

By successive differentiation of Eq.(|l6]) we can determine OA/dcj) = l/g{A) and d'^A/dcj)^. 

Evaluating the expressions ([T8| - pTD we obtain the results presented in Fig. 2, 3 in very 
good agreement with the numerical results, providing support for our approach [|15|. We 
also verified that the agreement remains (within a couple of percent) for other values of the 
interaction (e.g. U /t = 8). 

Comments: From our formulation of D{T) as a thermal average over level curvatures, 
we can try to understand the difference in behavior between integrable and nonintegrable 
systems. In nonintegrable cases the level repulsion prevents level crossings, i.e. crossings 
are statistically negligible for macroscopic systems. Then each level en{(p) fluctuates only 
on the scale Ae oc l/Af{E), where Af{E) is the many-body density of states. Therefore the 



curvature Dn averaged over or over different k vectors in the thermodynamic hmit should 
vanish. On the other hand, in integrable systems levels in general cross, so ffuctuations of 
e„(0) do not necessarily vanish for L -^ oo. Hence there is no restriction on the average 
D{T), except that D > 0. The difference between both cases is intimately related to level 
statistics. 

This connection between integrability and finite temperature charge transport born out 
of this model calculation, we can conjecture to hold true for other quantum (as well as 
classical) integrable systems. We can trace it to the existence of a macroscopic number of 
conservation laws. It is plausible that with respect to transport one should distinguish two 
types of models within the class of one-dimensional (ID) correlated systems where some 
solvable models [^] are available: 



a) In few (mostly solvable) ID models the current is a conserved quantity. This is generically 
the case for models without Umklapp scattering, e.g. the Luttinger model, ID Bose gas etc, 
but also for the U = oo Hubbard model [0. In these cases one expects at T > ideal 
conductance of the system characterized by p(T) = or D{T) > 0. 

b) Nontrivial answers are expected for models with Umklapp scattering, e.g. for the ID 
Hubbard model, the t — V model etc. If our conjecture is correct we expect integrable 
models as the Hubbard, t — V or supersymmetric t — J{J = 2t) models to behave as ideal 
conductors at finite temperatures. We can then argue that even ID nonintegrable models as 
the U — V model (with longer range interactions), as they are characterized by the integrable 
Luttinger liquid Hamiltonian at low energies, should behave as nearly ideal conductors at 
low T; the 6-peak then would broaden to a narrow Drude peak of weight D{T) |jl8[. 

We should stress that our study is based on the Kubo linear response theory Q] whose 
applicability in the context of Luttinger liquids has recently been debated [|19|. Further work 
is necessary on other systems to lend further support for these ideas. 
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FIGURES 
FIG. 1. D{T)/Dq as a function of th/t for U/t = 2, different temperatures T and system sizes 

L: (- • -) L = 6, (- • • -) L = 10, (- • • ■ -) L = 14. 

FIG. 2. a)lower and left axes: scaling of D{T = 2t)/Do with system size 1/L (circles) for 

tfi = t, U/t = 2. The indicates the Bethe ansatz result. The continuous line is a linear fit 

for L = 10, 14, 18 sites, b) upper and right axes: scaling of lnL'(T = 2t)/DQ with L^ (dots) for 

U/t = 2 and th = 0.5, 1.2, 2.0t 

FIG. 3. D{T)/Dq as a function of T for U/t = 2. Numerical results are indicated by 

(- • -) L = 6, (- • • -) L = 10, (- • • • -) L = 14 and th = 0.5t; by points for 

th = t and L = 18. The continuous line is the Bethe ansatz result and the dashed line D{T)/Dq 

for a free particle. 
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